PuDlic  reporting  burden  for 
maintaining  the  datanoede 
suggestions  for  reducing  this  burden,  to  wa»k« 
and  to  the  Office  of  Management  and  Budget  Papenvotk  Heduetton  rtvivs.^ 


PAGE 


Form  Approved 
0MB  No.  0704-018. 


Deluding  the  time  fut  icviewing  jnshocuvns  seai^ning  ex'simg  data  sources.  gath< 
s  regarding  this  burden  estimate  or  any  othei  aspect  ot  this  collection  of  information 
rpcratlons^dRcports.  t^lSJcIfersonDavisHiohwdy  Suite  1204  Arlington. VA  22202-4302. 
DC  20503  _ 


1  AGENCY  USE-ONLY  (leave  bfanA; 


2  REPORT  DATE 

August  1991 


3  REPOaTTVPS  AND  OATES  COVERED 

professional  paper 


4  TITLE  AND  SUSaUE 

AN  APPROACH  TO  THE  TARGET  DATA  ASSOCIATION  PROBLEM  USING 
SUBJECTIVE  AND  STAJ'ISTICAL  INFORMATION 


5  PUNOIHG  NUMBERS 


In-house  funding 


6  AUTHOR(S) 


I.  R.  Goodman 


7  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152-5000 


8  PERFORMING  ORGANIZAUON 
REPORT  NUMBER 


9  SPONSORINGyJ40NaORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152-5000 


10  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


■ELECriE^P 
SEP  3  5  1991 1 


n  SUPPLEMENTARY  NOTES 


12a  OlSTRIBUnON/AVAILABlLriY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


12b  DISTRIBUTION'! 
!?• 


13.  ABSTRACT  (Maximum  200  wo(d$) 


TR^T 


A  general  approacli  is  presented  to  tlie  data  association  problem  in  tracking  which  incorporates  both  linguistic- 
based  and  statistical  information,  extending  classical  results.  The  technique  continues  the  previous  effort  of  the  PACT 
(Possibilistic  Approach  to  Correlation  and  Tracking)  algorithm.  In  addition,  asymptotic  consistency  results  are  obtained 
relative  to  increasing  sample  sizes  and  model  refinements.  -  - -  -  - 


91-11482 


Published  in  the  Proceedings  of  the  1984  American  Control  Conference,  Volume  1  of  3,  1984. 


14  SUBJECT  TERMS 


-algorithm  / 


linguisticfbased  information 


probablistic  descriptions  ^ 


15  NUMBER  OF  PAGES 


10  PRICE  CODE 


17  SECURRYCLASSIFCATIOIJ 
OF  REPORT 

UNCLASSIFIED 


18  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


10,  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20  LIMITATION  OF  ABSTRACT 

SAME  AS  REPORT 


NSN 75400  S.28O-55C0 


Standard  form  203 


UNCLASSIFIED 


NSN754aOI.a80-5500 


Standard  form  238 


UNCLASSIFIED 


1984  AMERICAN  CONTROL  CONFERENCE 


i?  s 

>t.^  \ 

■  SMI 

*•  -j  .  L  . . '  t  .(  J 

UVi  *  .  •j 

_i' 

_  CU 

\ 

»•#  ^ 


I't 

-  ■*  .  J-  ••V, 

V-T—  ^  ^ 

I  I  .  :.  ;.  ■v‘-^-'T 

•  ■  C-  .i«  i  ^  ^ >  -.r  ■  ■  ■  ^' '  ■•■)>••{—••-  -.J  i  A;:  -i 


-y^  7- 


*■  ikr,  ■ 


AN  APPROACH  TO  THE  TARGET  DATA  ASSOCIATION  PROBLEM  USING 
SUBJECTIVE  AND  STATISTICAL  INFORMATION 


I.R.  Goodman 

Mm.  OCEAN  SYSTEMS  CENTER 
SurvoiUanoo  Systems  Dopartmont 
Code  7223 

San  Diego,  CAL.  92152 


Abstract 


A  general  approach  Is  presented  to  tho  data 
asscciation  problem  in  tracking  which  incorporates 
both  linguistic-based  and  statistical  information, 
extending  classical  results.  Tho  tochniquo  continues 
tho  previous  effort  of  tho  PACT  (Possibilistic 
Approach  to  Correlation  and  Tracking}  algorithm.  In 
addition,  asymptotic  consistency  results  are  obtained 
relative  to  increasing  sample  sizes  and  model  rofino- 
Ronts, 

T_.  Introduction 

"  ^  The  PACT  algorithm  has  been  documented  in  pre¬ 

vious  papers, ^l],[23.[^],qylt  Is  an  attempt ^to'  per¬ 
form  track-to-^r^iT'oT^eport-to-track  data  associ¬ 
ation  or  "correlation"  based  upon  varying  informa¬ 
tion  sources.  The  procedure  In  an  obvio,t/s  way  may 
be  extended  to  problans  of  medical  ^i^griosis,  fault 
determination,  and  classification  ,[!]}. ^he  entire  • 
technique  Is  based  upon  the  idea  thai  ^tributes 
may  be  represented  by  fuzzy  sets  and  tfiaYyUch  ob¬ 
jects  represent  a  weakening  of  probabil is'^c  des- 
criptions-and  hence  suitable  for  the  interpretation 
of  natural  language  informationCCSllW 'or  other  work). 

2.  Hybrid  Logical  Systems 
A  hybrid  logical  system  consists  of  the  tuple/® 
(^not* ^S'^or'*^^  whereJOnot  ^  collection  of  nega¬ 
tion  operators  9  :C0,1]-*[0,1]  nonincreasing  with 
<;j,(0)=l  )-C.  Other  properties  that  '5„'can  pos¬ 

sess  Include  involution  and  continuity  [5].  0.  is 
a  collection  of  conjunction  or  general  cartesian 
product  operators  0,: [0,1 ]‘*C0,1]  nondecreaslng^ 
continuous, bounded  above  by  min, and  such  that  for 
all  u,v:  <i,{0,v)=<>,{u,c)»0  ;  q.  (1  ,v)=v=<., (v.l ) ,  (1) 
as  well  as  higher  argument  operators  with* similar 
properties.  When  as  a  two  place  operator  is  sym¬ 
metric  and  associative,  then  it  is  called  a  t-norm, 
and  possesses  many  useful  properties,  including 
an  unambiguous  recursive  extension  to  multiple  argu¬ 
ments  and  cajionical  representations  [S  ],[6].  When 
n-arqument  0,  is  the  cumulative  distribution  funct¬ 
ion  (c.d.f.)*of  uniform  [0,1]  random  variables{r .v.) 
(not  necessarily  Independent, nor  need  be  symmet¬ 

ric  or  associative).  It  Is  called  S  copula  [6]. 

Copulas  and  t-norms  may  coincide  or  not. (See  [7].) 

A  similar  situation  holds  dually  for  C  ,  a  collect¬ 
ion  of  disjunction  or  general  cartesiafi'^sum  opera- 
Lors  ij*  :[0,lj^-[0,l],for  various  m,r.ondecroaslng 
continuous, bounded  below  by  max, and  e.g.,  for  m=2 
all  u.v:  Ojj(l  ,y)=<!  (u.l  )“1  ;  0  (0.v)-v=<i  (v,0).  (2) 

t-conorms  and  co-copulas  are  also  defined  analog¬ 
ously  [5]-,[6].  Finally,  It  Is  a  collection  of  fuzzy 
set  membership  or  possibility  functions  :X-C0,1  ] . 


corresponding  uniquely  to  attribute  or  fuzzy  subset 
A  of  space  X  having  domain  dom(A)»support(<>^)C X  ; 
for  various  A  and  X.  Such  functions  Include  ordinary 
probability  functions  v/hen  X  is  discrete,  ordinary 
edf  s  when  X?s'',unimodal  normable  continuous  or  not 
(nonnable= having- unit  value  somewhere),  etc.(See,e.o 
[9]  for  various  basic  properties  of  fuzzy  sets  and^ 
possibility  functions  as  well  as  many  applications.) 

The  collection  S  of  all  strings  produced  by  F 
consists  of  all  well-defined  combinations  of  operat¬ 
ors  from  F  evaluated  at  arbitrary  possibilitty  func¬ 
tions  chosen  from  U  and  possibly  further  evaluated 
at  functional  arguments.  A  typical  string  is  denoted 
comb( . - ,x ,..).( 3) 

The  following  principles  Involving  1 ingulstic/semsnt 
ic  evaTuatTons  are  assumed  ; 

Principl e  1  -  Parsing 

All  statements  in  natural  language  r^y  be  exp¬ 
ressed  in'aTonaal  language  in  terms  of  membershios' 
or  degrees  of  compatibility  of  elenents  attri¬ 
butes  and  the  use  of  operators  forming  formal 
strings  representing  (possible  different)  not, S. or. 
Symbolically,  a  typical  string  here  is 

s  “  comb(. .not. . .S. . ; .  .or. , ; . . (xeA) . .)  (a' 

=  ( ..  ,x,.  .)c(comb( .  .not. . ; .  .S. . ; .  .or . .)( .  .A. . )}, 

Principle  2  -  Abstraction 

The  truth  or  semantic  content  of  any  string  (or 
proposition)  as  in  eq.(A)  Is  tr(s),  as  In  eq.(3). 
If  the  system  Is  t.i  /'  ®.in: 'ional ,  then  (by  deflnlt- 

tr(s)=comb(..v„ -  .)(.., «^(x), . 

( ^ , 

Other  concepts  associated  with  natural  language  may 
be  described  and  evaluated.  Including:  conditioning, 
quantification,  interaction/independence  as  well  as 
dependence,  projections  and  marginal  attributes', 
tense,  mode,  case  and  prepositional  and  predicative 
relations;  and  functional  and  relational  transforms, 
Including  arithmetic  operations  and  linguistic  modi¬ 
fications.  (See  for  example,  [4]  or  [ 91.) Examples 
of  the  above  will  be  given  later  in  association  with 
tne  PACT  algoritJm  essentially  as  error  distribution 
and  inference  rules. 

One  example  of  a  hytirid  system  is  simple  proba¬ 
bility  logic  where  F»(neg;/j;t(^:L'^^^)  ,  where 

ncg(x)=1-x;\;>j  Is  a  fixed  n  (argument)copula ,  ^  Is 
the  OeMorgan°transform  of  and  thus  a  co-cop8la. 
and  is  the  collection  of  all  cdf's  over  R  . 

Then  by  Sklar’s  Theorem  [7]  for  any  F.  ,..F 
■?,(?,  ...,F  )  (In  compositional  form)  Is  a  s?.mple 
string  repPesenting  a  legitimate  cdf.  By  replacing 
Fj's  by  1,  marginal  cdf's  or  cdf's  of  number  less 


than  n  can  ie  fonned.  More  generally ,ron-cdf  string 
coab(neg;<»j,i:(^)(F^  ’  •  ‘ 

•  . , F^)(co(nb{C,X,t )({-". ^1  ]*  •  •  . 

which  may  be  extended  uniquely  abcve  the  ray  level 
to  arbitrary  Borel  sets,  but  for  the  latter  situa¬ 
tion,  F  is  no  longer  truth  functional  in  general 
(a  characteristic  of  jointness  of  events  in  prob¬ 
ability-  see,  e.g.  Rescher  [10]). 

More  generally,  a  Lebesguc-Stieltjes  logic  may 
be  developed,  where  cdf's  are  replaced  by  functions 
.of  bounded  variation  and  probability  measures  by 
signed  measures  at  the  Borel  set  level  [?  ].  Of 
course,  in  the  probability  logic  example, ^he  single 
if,  could  be  replaced  by  a  collection  of  copulas  and 
related  functions  (noting  that  in  general  copulas 
composed  with  multiple  copulas  are  not  copulas, 
although  they  do  rejnain  conjunction  operators). 

A  single  natural  language  concept  may  lead  to 
a  variety  of  possible  semantic  evaluations.  One 
contributing  factor  to  this  ambiguity  lies  in  the 
nonuniqueness  of  extending  an  ordinary  function 
frX^xX^^Y  to  f:F(X^)KF(X2hF(Y),  the  lifting  of  f. 

where  r (X)  denotes  an  appropriate  class  of  attri¬ 
butes  of  X  (  or  all  attributes  on  X).Two  candidates 
for  ?  are,  using  Principles  1  and  2  and  the  unique 
representation  of  attributes  through  their  member¬ 
ship  functions  for  a  simple  logical  system; 

o>(y)=*r(ycf(A,B))=tH(3(x,2)£A*B)i^“f(x,i5^ 

(all  (x,2)ef‘'(y)) 

and  .  . 

'>5;^,B)(y)=tr{f"  (y)cAxB)“<!.j(.jA»^B^^^" 

where  in  the  first  case,  AcF(X, ),Bc?(X-) ,  and  in 
the  second  case,  AcF(P(X-,)),  Bef(P(Xp))T  related  to 
the  measure  theory  concept  of  transformation  of  mea¬ 
sure. 

In  any  case,  once  the  functional  extension 
problem  is  determined,  other  ambiguities  arise:  Let 
SP  be  the  class  of  all  multiple  and  single  argument 
Boolean  truth  polynomials,  i.e., 

SP 'U({0,l}^°'’n'  (9) 

n=l 

Thus  each  function  f  in  SP  may  be  extended  to  f: 
F[0,l3y..xf’[o,l Vp[0,l].  Call  this  class  S?>.  flow, 
f=coaib(?^;<(j;<.^)  restricted  to  {0,1)”  for  suitable 

choices  of  n  yields  a  classical  logical  connect¬ 
ive  ,  and  in  turn  the  extension  f  thus  provides 
an  alternative  natural  evaluation  to  ea.(5) 

tr(s)x(<f^(  ^,^A)..)^’°>*’>f{..(A,CA)..)^’^^- 

in  general,  distinct  in  form  from  that  in  eq.(5). 

Another  problem  is  the  choice  of  the  operators 
for  F  to  .'node!  a  given  si..uation.  Restricting  these 
operators  to  be  ,  e.g.,  t-norms  or  copulas  may  be 
confining,  despite  the  appeal.  This  is  closely  re- 
-lated  to  Hilbert's  13Lh  problem,  and  more  generally, 
to  the  problem  of  representing  multiple  argument 
functions. by  superpositions  of  lower  argument  func¬ 
tions  Dl3,v/hen  associative  copulas  or  t-.nonr.s  are 
employed,  because  of  their  canonical  representations, 
Gn  the  positive  side,  it  can  be  shewn  that  SprechcHj 
refinor.ent  of  Kolmogorov's  original  form  for  an  arb¬ 
itrary  continuous  function  of  multiple  arguments  may 


be  formally  couched  as  the  strino^for  constants  Cj, 

wnere’ Vndsum  '.s  ordinary  sum  bounded  unity, 
a  t-conorm,g  is  monotone  decreas.ng  ,  and  ■!>,  is  a 
strict  (  and  hence  Archimedean)  t-norm  having  can¬ 
onical  representation  , 

,..X|^)=h"‘(h(x.|)+..+h(x^)),  (12) 

where  hrfO.l  is  continuous  decreasing  with  h(0) 
o+®  and  h(l)=0.  Also,  the  X-intensification  or  mod¬ 
ification  of  A  through  is  given  by 

4^^Ua^x))^  h-''(xh(4^(x))),  (13) 

noting  that  these  intensifications  obey  the  exponen- 

'*•> .  (i«) 

for  all  nonnegative  real  X,  extending  the  integer- 
valued  repetition  of  operator  <»,.  A  reasonable  sub¬ 
stitution  for  the  above  form  is^  the  use  of  ordinary 
exponentiation  or  translation  of  arguments.(See[’3  3J 

2-  Combination  of  Evidence 

Let  F  be  a  given  hybrid  logical  system  to  be 
used  for  the  relevant  modeling.  Let  e  c  X  be  an  un¬ 
known  parameter  vector  to  be  estimated  or  decided 
upon.  Suppose  that  :XwYj*[0,l],  ,*l,..,m,  are  m 

possibility  functions'^describing  s  as  well  as  intro¬ 
ducing  nuisance  parameters  Zj  cYj  ,  j=>l,,.,ra. 

Then  ,  if  for  example  a  single  m-copula  i.  is  chosen 
for  combination,  one  has  the  joint  non-interactive 
function  -  ~ 

«^(e,Z)»<<j(9^  .  (15) 

0^(e,....8^)  ;  Z^(Z,...,Z^)  .  (16) 

This  representation  of  the  joint  information  and 
all  nuisance  values  may  be  intercreted  in  terms  of 
either  ordinary  operations  upon  nested  level  sets 
or  equivalently  upon  cartesian  products  ?  in  partic¬ 
ular,  on  nested  random  sets.  (Indeed,  a  strong  iso¬ 
morphic-like  relation  holds  between  many  common  set 
operations  on  nested(randcm)sets  and  corresponding 
possibility  functions,  including  intersections, 
unions,  set  differences,  projections,  and  functional 
transforms.  See  [6]  and  [12.].)  Thus, for  all  e,Z  j 


,j(6.Z).Pr(to  Z,  )=S  [A,  S5. 
where  S,,  (A,)2^  [U,,lJ  is  the  unique  nested  random 
“j  '^j  ^ 

subset  of  XxYj  which  is  one  point  coverage  equival¬ 
ent  to  <f.  ,  where  U,  ,..U  are.marginally,rv's  dis- 
tributecj  J  uniformly  over  [0,1]  with  joint  cdf  o,- 
(See  [6], [7]  for  other  properties  and  results.; 

Next,  consider  the  concept  of  conditional  pos¬ 
sibility  functions:  Let  X  and  V  be  given  spaces 
with  attribute  D  c  f(Xxy)  dominated  by  attribute 
C  c  FiV)  i.e., 

sup  *n(x,y)<  ^Ay)  \  all  ytV  .  (1^ 

(xcX)  “  ^  ^  ^ 

Then  it  foUuws  that  there  exists  (uniquely,  if, eg. 

is  a  strict  t-norm)  conditional  attributes(Blv) 
satisfying  (for  given  operator  <'j)for'  all  xcX.yc/ 

<>g(^  >y)''l'j{^{jj|^|^)  ■  (19) 

in  particular,  let  the  diagonal  possibility  doscrip* 


5R8 


Urtn  of  0  be  \  assuming  here  a  simple  system^ 

(all  Z) 

Let  X  -X  and  /  be  the  space  of  all  possible  descr.ip- 
Ll'bns  (  or  ’’possible  worlds")  of  8.  Hence, d{A)£y 
and  the  universal  Inclusion  description  Ccf(y)  with 


d  (A ) ) = tKQ  d )  ( 0  ed  (A  )))=o 
^  .'f- 


i^d(A) 


(0)). 


(21) 


Since  the  compound  attrlbote  representing  o  and  d(A) 
has  the  same  evaluation.  It  follows  that 

«d(A)(0)“<'s(«(®l‘l(A)).<>p(d(A)))  (22) 

determines  the  final  combination  of  evidence  func¬ 
tion  i!i(old(A))  as  a  function  of  OcX. 

As  a  check  on  this  procedure  reducing  to  the 
classical  combination  of  statistical  evidence,  if 

c>.=prod,<>  *bndsuo,Y.  Is  vacuous,  and  X=ll''suitably 
&  or  J 

Hlscret-’zed  where  each  Is  approximated  by 

“^A  (Oj)‘“^(°j"“j*Aj)*A(0j)  .  (23) 

V  being  the  Standard  probability  density  function 
(p.d.f.)  for  the  gausslan  distribution  1^0,1),  then 

a(oJd(A))=l'(o-u,A)*  a(0)  ,  (2A) 

a-^aJ  -Uj  ;  A^(’SA'’)‘^  .  (25) 

the  same^firoally  as  th4“kuE  (best  linear  unbiased 
estimator)  of  u,and  by  a  fiducial  argument,  of  e. 

£.  Application  ^  Data  Association 

Suppose  now  that  attributes  describ¬ 

ing  nuisance  parameters  concerned  with  e  (  but 
not  e  directly)  are  available  In  the  form  of  com¬ 
pound  attributes  represented  by  error  functions 
and  by  Inference  rules.  A,  could  represent  geolo¬ 
cation,  A»  a  sensor  systen  state,  A^,  another  state, 
and  for  mM,  A,  could  represent  vlsOal  descriptions, 
where  o  represents  the  level  of  correlation  between 
a  pair  of  track  histories,  say  1  and  J.  Thus 


described  separately  by  error  functions 

fori*i,j  and  k*l,..,n  ,  where  are  observed  or 
predicted  outcomes  of  data  associated  with  A.  -  In 
fact,  1^  e  do(n(Aj.)  -  and  c  dom(Aj,)  can  be 

arbitrary.  Denote  ,■ 

Z=(Z^...,Z^)  ;  Z^=(zn),Z^>l))  ,k=1...,n  (26) 

with  similar  notation  for  Z  and 

In  addition,  Z.and  8  are  described  jolntlyCor 
perhaps,  more  accurately,  conditionally)  by  r  Infer¬ 
ence  rules  R. ,  t=l,..r,  where 

<>{^(0|  2)“  tr((Z  c  G^)  iff(o  £  Corr^  )) 

‘  ■  (01,  '  (27) 

t  Y- 

where  Implication  operator  or  <»  ,  corres¬ 

ponding  to. "If-then"  and"1ff"'^,  respectively.  For 

(28) 
r) 


example,  one  can  define  for  all  u,vc[0,l],, 

«♦«(  u .  V ) -« ,  ('>^(u ,  V ) ,  ,u ) ) ;  Vj,(  u .  V ) -^^  ( 0^  { u ) .  V ) 

where  antecedent.  ,  /  » 

if.  (Z)"tr(  £  (Z/^  t  zM  match  to  degree  6.  ^)] 

(28) 


fkcJj.)  ®t,k^ 

J.9{l,..,m)  Is  index  set  of  attributes  operates 


on,  M.  is  a  matchIngSet  for  attribute  A.  so  that 
^^k^  degree  to  which  Zj,-  i  and  Zj,lJ' 

are  compatible, with  Intensification  (or  extensifl- 
catlon)  dctenalned  by  iip  for  some  j^aO  , either 

®t,k 

In  the  form  as  given  in  eq.(13)  or  (as  previously 
mentioned)  as  exponential  or  translational  argument 
form.CorrY,  similarly  represents  y-deoree  correla¬ 
tion.  A  tabne  can  be  established  relating  linguist¬ 
ic  intensifications  with  numerical  ones -e^j,,Yt. 

From  now  on,simpl1fying  notation  for  all  rele¬ 
vant  error  and  matching  functions  and  Inference 
rules  win  be  employed,  basic  problem  here  is 
to  apply  the  combination  of  evidence  procedure  to 
the  data  association  problem  for  an  appropriately 
chosen  hybrid  logical  system. First  note  that: 

A.  Attributes  may  be  grouped  roughly  Into  two 
parts:  (A-j .  ,A^, },  statistical' attributes  ; 

(A  A  },  subjective  attributes . 

For  each  statistical  attribute  such  as  geoloca¬ 
tion  or  sensor  state  parameter,  ^Pj,  arises  typical¬ 
ly  as  the  discretization  of  pdf  fj^,  say,  obtained 
by  various  geometric  and  physical  considerations. 
For  each  subjective  attribute  A.  such  as  Identifi¬ 
cation  or  visual  description  of*‘some  type,  Pj^  Is 
obtained  from  a  panel  of  experts  and  represents 
possibility  values  conditioned  on  potential  observ’- 
ed  data  and  In  general  will  be  approximately  a 
symmetric  binary  function,  so  that  a  reasonable 
choice  for  Hj,  1n_any  is  F'j,(Zy,)'P]j(2^  M^^^)). 

B.  On  the  other  hand,  for  any  statistical  attri¬ 
bute  A.  ,  In  general  M.  is  not  the  same  as  Pj.,  since 
statistical  hypotheses  testing  (Involving  P^.)  pre¬ 
sents  a  natural  candidate  for  H.  :  the  significance 
level  function  for  one-dlmenslorial  test  statistic 

say,  where  nj.'‘nj,(4j,)  for  any  outcome  of  rv  Zj,, 

noting  anatn  that  P.  is  a  discretization  of  f.  the 
pdf  oflZ^Z^.  Thus, for  the  hypotheses 

HqII.J  match  vs.  :  1 , j  don’t  match  ,  one  obtains 


«l(2,)-Pr(,^(5,)>,^(Ij)|K  ).l-T,(,^(Z^)).  (29) 

where  Tj,  is  the  cdf  oF(n;^(z^)!HQ). 

C.  Inference  rules  are  generally  selected  from  a 
panel  of  experts  and  arc  thus  uniquely  determined 
by  -once  a  choice  of  F  Is  made-  0*  and  the  tuple 
((Bt  lY^)  .  although  nodlflcStions  to  this 

may  occur  If  negative  matching  or  negative  correla¬ 
tion  relations  are  Introduced. 

The  following  three  theorems  will  be  useful  In 
developing  principles  for  choosing  appropriate 
classes  of  operators  for  F. 


Theorem  1 . 

Let  BeX  with  tF(X),  Vj  a  conjunction 

operator,  agd  form  for  any  possible  value  x  of  e,. 

pj  0  (*)“p5(^g 
k-1  2  ^  ^ 

Then:  ’ 


If  Is  either  an  Archimedean  t-norm  or  min 
6 


(xlxcX  fi  11m  (x)»ll-  /5 


k— +-  k 

then,  unless  •  min,  «q(x)  •  0,  for  all  xcX. 


(31) 
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(43) 


fiJ)  The  a-1eve1  description  of  Q  by  level  sets  equ¬ 
ivalent  to  Bv,02.--  Is  , 

g{a)=  n  [a.l])cX  ;(^o<l  .  (32) 

l<k<-s»  “k 
Then  i 

«n‘'[a.l]Cg(a)  .  all  Osft<l  ,  .  (33) 

with  proper  subset  relation  holding  in  general, and 
with  equality  holding  iff  <{•,  =  min.  i.e.,  the  only 
case  where  no  information  is  lost  concerning  e  is 
min. 

Theorem  2. 

LetHCcHR  )  vfith  the  discretization  C  c 
f(D  )  of  C  defined  by  letting  Dp  be  the  domain^of 

C  formed  by  making  a  discretization  of  r"  with 
mBsh  fineness  strictly  decreasing  as  p  grows, etc., 
and  where  ^  ^34) 

t  i-  P 

Then,  ^ 

(i)  If  I}’  is  an  Archimedean  t-conorm, 

°  lira  i  («  (x)}  =  1  .  (35) 

^  (xcD„)P 

(ii)  If  =  max  , 

°  lim  (x))  il  ,  (36) 

0 

with,  in  general,  strict  inequality  holding. 

Theorem  3.  (Soe  also  Cl 33.) 

(i )  Let  •*■1^  .  k=l,..,a, for  affixed  finite, 

be  bounded  pdf's. 

(ii)  Let  '{'.•jbe  a  conjunctioa Define  the  p  modffied 
discratizatiin  fj,  ^  of  fj.  as 

'^^k  '  °k,p' 

where  Or  is  the  corresponding  discrete  domain 

K ,  P 

and  a  is  the  mash  function  so  that 

and  finally,  let  for  any  x=’(x^ , .  .,x^)  c  O^pj  - 

^(p)^^^''^?/^l  .p^^l  ^  •  ’^a  ,p^^a^ 

(iii)  Suppose  also  that 

’*’jl^“l'"'“a^'^^“r'^^'a^u,“-.=u  =0 
exists  finitely.  ^ 

(iv)  Suppose  that  <{>,-  is  a  conjunction  such  that 
there  is  6,  0<6<1,  sucn  that  for  all  ucCO.l], 
3i;fj2(u,v)/3v  and  3^^^(v,v)/5^y  are  bounded  for  all 

vc[0,6). 

(v)  4_:X><1  ®  -»  [0,1]  is  a  given  continuous  function 
v/here  *X  is  the  space  that  contains  0  and  where 

m.  ^  m,+. .+m,  .  (41 ) 

O  I  a 

(yi)  It  is  assumed  that  <>  is  an  Archimedean  t-co¬ 
norm  with  canonical  generator  h.'h  is  continuous  with 

$^(u,v)°l-h"^(min(h(T-u)+h(l-v) ,h(0)))  ,  (42) 

for  all  u,ve[0,l3,  snd  h:[0,l]-*(R^  is  nonincreas¬ 
ing  with  h(l)'‘0.(See  [5]  for  further  properties.) 

It  is  further  assumed  that  there  is  a  <,  0<«:<l,such 
that  dn(v)/dv  and  d^h(v)/d^v  are  bounded  functions 
of  vc[ir,13. 

Then,  for  all  OcX  , 


lim{  (<>j2('>i-fo«’')-^(p)(*5))  ) 

-  l-h"^ (min(T{8) .h(0)))  , 
where  3 

^  xclf°  ^  k-1-''- 

fe(u)^-{dh(v)/dv)^,^  •(3«J.j2(u,v)/Dv)y^g  ,  (45) 

for  all  ucC0,l3. 

Remark. 

Frank's  family  of  modular  t-norms  (see  [b  J  or 
[6  3  for  various  properties )and  t-conorms  can  be 
shown  to  satisfy  the  conditions  of  Theorem  3  for 
any  ,  <l>^2  ^nd  chosen  from  the  Archimedean 

part  of  the  family. 

Theorems  1,2,3  lead  to  the  following; 

Principle  3 

If  <.g  :X-C0,13  ,  k»l,2...  is  a  sequence  of  des¬ 
criptions  Sf  0gX  with  either  high  redundancy  or  ^ 
much  irregularity  occurs  in  the  sequencr 
lim  (0)  -1  does  not  appear  viable,  t« 

{>*4«  “k 

for  combining  this  sequence  i^j^nain.  In  addition, 
min  allows  for  an  uncountable^nicnber  of  descriptions 
of  0. 

Principle  ^ 

If  as  above  are  now  relatively  non-inter-* 

active/non-overlapping,  then  <!»,  and  4  as  Archimed¬ 
ean  t-norms  and  t-conorms,  respectively,  may  be 
reasonable. 

Principle  £  .. 

If  C  is  the  discretization  of  C  as  in 
Theorem  2?  then  for  the  model, choose  e-  =max  so  that 
Hm  (  ifohr  becomes  non-trivial 

p-f+»  J^cDp  p 

Thus,  returning  to  the  basic  problem.  Principle 
4- Implies  that  if  A, , .  .,A^,8re  all  relatively  non¬ 
overlapping,  then  the  total  error  effect  may  be  ob¬ 
tained  by  simply  choosing  an  appropriate  conjunction 
Y<hich  may  well  be  an  Archimedean  t-norm  such  as 

from  Frank’s  family.  This  yields  the  overall  stat¬ 
istical  error  effect  . 

P(Z’!2’)i  i-,  ,(Pt(ZJZJ)  .  (46) 


'  ’'J.r"k'‘-k'‘-k^ 

dsksm') 


for  all  k,  k^l  ,. . ,ra. 

Similarly,  for  .  tb.e  overall  sub¬ 

jective  error  effect  is 

P(Z-|Z")^  •  '•'8) 

('likSml 

Y/Ith  possibly  the  two  conjunctions  being  chosen  the 

same, and  further, as  an  Archimedean  t-norm. 

On  the  other  hand,  since  inference  rules  typic¬ 
ally  operate  on- combinations  of  common  or  similar 
attributes  -  by  intensification  modifications-  Prin- 
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ciple  3  may  be  invoiced  to  yield  the  overall  infer¬ 
ence  rule  effect  as 

RC0lZ)=niin  (R/'oiZ))  .  (4S) 

l<t<r 

Again  using  Principle  4,  it  ncay  be  eppropiMate 
‘to  choose  Arch.siedean  t-norm  bo  obtain  the 

joint  posterior  function 

^{O.ZIZ.R.P)^  «^3(p(zi!).r(o1z)) 

o  <.^3(P(Z'l?’).«i.(0.Z))  ,  (50) 

Khere  overall  error  effect  is 

P(Z|Z)^<.j3(P(Z‘|Z’).P(Z“|Z”))  (SI) 

and, as  in  Theorem  3, 

^j.(0,Z)^«j3{P(Z“|2-).R(0|Z)).  (52) 

Note  that  a  sufficiency  condition  holds  here 

«(8lZ,Z.R.P)-«(0lZ.R)=R(6lZ)  (S3) 

as  well  as  -fct-  „  o  j 

«(Zl?.R.P)=«(Z!2.P)=P(Z'Z)^  (54) 

which  also  when  substituted  into  (50)  yields  a  form 

which  may  be  independently  derived  by  use  of  con¬ 
ditioning  of  possibility  functions. 

Hext,  Principle  5  may  be  used  to  obtain  the 

marginal  posterior  function  for  6  ,  before  normali¬ 
zation,  using  the  universal  inclusion  description: 

-  (  i'_(4(6,Z|Z,R,?))),  (55) 

'  fell 

for  some  choice  of  Archimedean  t-conorm  .followed 

‘in  turn  by  again  Invoking  Principle  5, 

4n(d(A))»  sup  {<S(ei2,R,P))  ,  (56) 

O<0<1 

which  finally  yields  the  normalized  posterior 

f(e|d(A))  (dependent  upon  2,R,P)  ,  by  use  of  (22). 

Note  that  if  P.  (Zj,|?j.)“<5-  2  (Kronecker  delta), 

^  ^  ^ 

then  the  evaluation  of  all  R^  involving  that  attrl- 

bute  A.  ,in  effect  have  M.  (Z. )  replaced  by  H.  (Z. ) 
in  the^antecedent  and  the  -term  Pi,(Zj^Iz. )  is^dropped 
from  the  calculations  for  the(norm3rized)  posterior. 
In  the  above  computations,  for  simplicity,  the 
for  Z=l,2,3  could  be  set  equal,  if  other  conditiSns 
do  not  arise. 

5. Asymptotic  Behavior 

Suppose  that  the  combination  of  evidence  prob- 

lan  relative  to  data  association  is  carried  out  as 

outlined  in  the  previous  section,  where  P.  is  the 
p'-"  discretization  (modified)  of  pdf  f,^,  for  kal,,.. 

m' ,  and  all  corresponding  conditions  for  parts  (i)- 
•  (iv)  and  (v1)  for  Theorem  3  hold  (a=Zm').  Also, 
note  that  part  (v)  of  Theorem  3  is  satisfied  by  use 
of  eq.(52)  and  that  any  given  inference  rule  R.  may 
use  Implication  operator  -►  being  either  ♦  or^>  . 

(A)  Under  the  above  general  conditions,  for  all 
OcCO.l]  . 

Hm  («(o|i.R,P)) “  l-h‘^(mln(Tlo).h(0)))  (57) 


as  in  eq.(43),  but  changed  to  account  for  max,  where 
t{6)'=/i,„.-  max  (  E,.(H«p{0.(V.Z"))))  ).  (58) 

where  I  in  eq.(52)  is  replaced  by  (V.Z")  ,  V  being 
a  r.v.  representing  formally  Z’  and  having  pdf  f: 

f(Z'|?-)-  f’(f.(z^’^|z^’^)-fj^(z^-^^|z[-^^)). 

3c  k  Jc  fc  k  K 

(B)  If  also  both  statistical  and  possibilistic  con¬ 
sistency  holds, l.e., 

-  <5(zj['^^-?[‘^^)  (dirac  delta){60) 

for  k»l,..,m‘  ,  and 

p,(z,W|5(«)  -  6z^^)  (Kronecker  dlt3)(61) 

for  k=m'+l,..,m  ,  and  if  k  is  such  that  k  is  non- 
decreasing  with  fe(0)=0  (  as  is  true  for  the  members 
of  Frank's  family-see  [?)  ),  then 

and  hence 

11m  (^(e|Z.R,P))  <  l-h"’(A.  .-fed)).  (63) 

p-,-+«  - 

(C)  Suppose  the  conditions  in  (A)  and  (B)  obtain. 

Then  o 

R(0i5)<  min  (i^„(l-6.(Z),9Tt)) 
litsr  °  ^ 

<  min  (  card(jj.)-^^j^(Z.^)®t.k)+  e^b) 

S-lff  (60 

(0)  Thus,  if  also  for  k=l,.,,m  , 

constant,  (65) 

then  o 

R(o|Z)<  min  (card(j  )(l-(l-<a)®t)+  e^b) 
lst<r 


<  m-(l-(l-a)®o)  +  6^0 


where 


Bji  mx(«  )  ;  9  i  4  mx<TO  .  (67) 

Hence,  for  all  ec[0,l]  , 

lim  (  «(0|Z,R.P))^-  1-h”^{ft2p.-fc(co)).  (68) 


Thus,  if  also  H,  (Z|.)  -  1  uniformly,  for  k”l,.. 
m,  is  bounded  positively  below,  and  Increases 


without  bound,  then  it  follows  from  (66)  and  (68) 
that  fer  all 

lim  (  «(e|Z.R.P))  «1-h"‘(0)=0.  (69) 

Finally,  since  in  general  i>p(d(A))>0,  (22) 
implies  that  the  final  normalized  posterior  is 

lim  (  «(0|d(A)))  °  ^;all  OcCO,l],(7o) 

p-»+—  ’ 

This  may  be  interpreted  as  describing  the  correla¬ 
tion  condition  'extremely  likely  correlated", since 
this  can  be  considered,  e.g.,  a  limiting  case  of 
the  function  «^(o)=o‘^  ,  for  all  o,  as  0“+®,  where 

all  intensification  Is  by  exponentiation  applied  to 
a  standard  membership  function- in  this  case,  the  1- 


dentity  function  over  [0,1]  representing  the  attri¬ 
bute  "correlation  level"  with  domain  being  probabil¬ 
ity  of  correlation  or  possibility  of  correlation. 

(E)  On  the  other  hand,  with  the  conditions  in  (A) 
and  (B)  holding,  the  lower  bound  for  R  is 

R(0if)>  inin(<-,{6  (1-G  (2).0^t),6  (l-oift.G  (2)))) 
^lSt<r)^  °  ^  °  ^ 

>  min  (  I-GJD-O^t) 

<1  st<r)  „  „ 

>  rain  (1-  min(H.  (Z. 

~(Ut<r)  (jccJ.)  ^ 

«P„.  (71) 

(F)  Thus,  if  also  for  k=l,..,o 

ti)iM^(?v,}>0  ;  0<fi)<l  constant  ,  (72) 

then  „  .  . 

R(o|Z)>  min  (l-u”t  -e^t) 

(l<t<r) 

a  l-u®o'-07o' 


where 


8  roinfS.  ,)  ',8  min(6^'5',Y„'-  minCrJ 

t. 


Hence,  for  all  o  c  [0,1]  , 

lira  (  «{0lZ.R.P))>  l-h'\/u„.-fe{o J).  (75) 

Thus,  If  also  •*  0  uniformly,  for  k=l,.., 

ra,3  ‘is  bounded  positively  belcw,  and  y'  Increases 
witnout  bound,  then  it  follows  from  (73°  and  (*75) 
that  for  all  0<9<1 , 

lira  (  «{0iZ,R.P))  >  l-h"’(/i.  .•fc(l))  .  (76) 

p-H-co 

But,  eq.(63)  shows  the  other  direction,  and  hence 
Hni(  ^(eIZ.R,P))  =  .-fed))  >  0  (77) 

in  general  (as  is  true, e. a.,  for  Frank's  farailyt^]). 
Also  note  that  for  8=1  , 

R(l|f)=nin  (G,(Z))  <  u^o  ,  (78) 


where 


n  it  <r.  V 
\R^  uses<>) 


Sq  -  n1n(S  j^) 


Acdt- 

\Utfr; 


R^  uses 


Thus,  if  there  is  an  R,  using  4-^,  it  follows 

that  if  B"  is  bounded  positively  from  below,  then 
the  above  conditions  imply  that 

lira  (dOliR.P))  =l-h-^0}=0  ,  (80) 

p-»+“ 

and  hence  from  (22)  it  follows  that  for  all  o  c 

[0.1], 

lira  (  4(o|d(A)))  -  1  -  <5  ,  ,  (81) 

p^4«  ' 

which, analoguous  to  the  result  ln(0),  indicates 
linguistically  that  the  condition  "not  extremely 
correlated"  holds,  corresponding  to  the  limiting 
case  for  possibility  function  ^  (o)=1.0R  as  a  func¬ 
tion  of  0,  as  ® 

In  summary,  the  above  results  show  that  suf¬ 
ficiently  good  Information  and  matching  of  data 
km'"  ^  leads  to  asymptotic  consistent  results, 
while  relatively  poor  matching  leads  to  the  unde¬ 


sirable  results  of  (F). 

6.Concludinq  Remarks 

This  paper  has  shown  how  a  combination  of  evi-  / 
dence  procedure  may  be  established  with  a  general 
application  to  the  data  association  problem.  Some  -it 
asymptotic  consistency  properties  of  the  procedure 
were  established.  Imp! c-mentat1xni  of  the  procedure 
Is  necessarily deaply  involved  with  the  use  of  Kalman 
filters  for  updating  geolocation  state  vectors,  as 
well  as  other  real-world  procedures.  Including  gT 

sensor  system  models  and  schemes  for  the  clustering  W 
of  ordered  pairs  based  upon  their  level  of  correla-  j/- 
tion,  often  presented  In  the  form  of  a 'correlation  !|j> 
matrix”.  The  latter  depends  upon  a  single  figure-of- 
merit  representing  the  average  correlation  level  (or  i§r 
probability)  between  any  two  track  histories  i,j.  !p 
One  suggestion  for  this  has  been  to  use,  in  effect.  , 
the  two  stage  r.v.  (H|Sj.(C))  ,  where  S..(C)  is  the  -y 
naturally  correspond! ng^random  set(d1scussed  before  -y 
-see  section  3)  for  output  4  ’■<’(*  ld(A)) ,  a  function  I' 
of  6e(0,l],  as  obtained  abovi,  and  where  H  in  its  }V 
conditional  form  is  a  r.v.  uniformly  distributed  y 
over  the  outcome  S|.(C)  [14].  In  particular.  If  ;i- 

<>j.  Is  uniraodal,  then  the  measure  of  mean  is  J:.' 


E(‘^lSy{C))=E(4/2)(4^ 

f*  .» 


■’(u)  ) 

eft  branch  right  branch 


=1/2  (  }  xd^Jx)  +  !r  xd(l-4^(x))) 

x=0  ^x=u  .  j-; 

’^vi’Av)  +  (l/2)(  /4p(x)dx  - /4p(x)dx  (82)  \' 

-  ~  x=u  x=0'-  ^  !-• 

where  u  Is  the  mode.  An  open  question  connected  :[ 
with  this  value  is  the  determination  of  a  related  !-/ 
figure-of-merit  which  is  invariant  with  respect  to 
any  particular  random  set, one  point  coverage  equi¬ 
valent  to  .  ■ 
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